Inner Accretion Disk Edges in a Kerr-Like Spacetime 
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According to the no-hair theorem, astrophysical black holes are uniquely described by the Kerr 
metric. In order to test this theorem with observations in either the electromagnetic or gravitational- 
wave spectra, several Kerr-like spacetimes have been constructed which describe potential deviations 
from the Kerr spacetime in parametric form. For electromagnetic tests of the no-hair theorem, such 
metrics allow for the proper modeling of the accretions flows around candidate black holes and the 
radiation emitted from them. In many of these models, the location of the inner edge of the accretion 
disk is of special importance. This inner edge is often taken to coincide with the innermost stable 
circular orbit (ISCO), which can serve as a direct probe of the spin and the deviation from the Kerr 
metric. In certain cases, however, an ISCO does not exist, and the inner edge of an accretion disk 
is instead determined by an instability against small perturbations in the direction vertical to the 
disk. In this paper, I analyze the properties of accretion disks in the Kerr-like metric proposed by 
Johannsen & Psaltis (2011), whose inner edges are located at the radii where this vertical instability 
occurs. I derive expressions of the energy and axial angular momentum of disk particles that move 
on circular equatorial orbits and calculate the locations of the inner disk edges. As a possible 
observable of such accretion disks, I simulate proflles of relativistically broadened iron lines and 
show that they depend significantly on the values of the spin and the deviation parameter. 

PACS numbers: 04.50.Kd,04.70.-s 



I. INTRODUCTION 

According to the no-hair theorem, black holes are 
uniquely characterized by their masses M and spins J 
and are described by the Kerr metric (e.g., [l|)- For 
this reason, all astrophysical black holes are expected to 
be Kerr black holes. Evidence for the validity of this 
assumption, however, has so far been inferred only in- 
directly (e.g., |2|). It is, therefore, desirable to test the 
no-hair theorem through observations. 

Tests of the precise strong-field nature of black hole 
candidates are often constructed in a model-independent, 
phenomenological approach, which does not assume any 
particular theory of gravity, and the underlying theory is 
generally unknown. Instead, these tests encompass large 
classes of modified theories of gravity and insight in the 
fundamental theory is hoped to be gained through ob- 
servations. It is then assumed that particles move along 
geodesies, which allows for the calculation and prediction 
of possible observable signatures of the theory. The as- 
sumption of geodesic particle motion is consistent with 
current experiments (see Ref. J3l]) and holds in general 
relativity and in at least all modified theories of grav- 
ity with known black hole solutions that differ from the 
Kerr metric [4] . Several such strong- field tests have been 
suggested, which can be performed with observations in 
both the gravitational-wave |5l-[l5| and the electromag- 
netic spectra |16l - [25j . 

These strong-field tests rely on parameteric spacetimes 
that deviate either slightly or severely from the Kerr met- 
ric, and several such Kerr-like metrics have been devel- 
oped to date (e.g., @, S IH El HI HI). See Ref. [H 
for a detailed discussion of the strong-field properties of 



these metrics. Since general relativity is still practically 
untested in the strong- field regime [29|, deviations from 
the Kerr metric, should they exist, could be either small 
or large as long as they are consistent with the limits 
imposed by current weak-field tests of general relativity 
(see Ref. Q). 

If a non-zero deviation from the Kerr metric is detected 
by either a gravitational- wave or an electromagnetic mea- 
surement, there are two possible interpretations. Within 
general relativity, the only possibility is that the com- 
pact object is not a black hole but instead either a stable 
stellar configuration or an exotic object [^, |30| • Alterna- 
tively, black holes are different from the expected Kerr 
black holes (see, e.g., Refs. |3ll - [33| ). and general relativ- 
ity is modified in the strong-field regime. 

Other tests of the no-hair theorem exist in both the 
gravitational-wave and electromagnetic spectra. After a 
merger with an object of comparable mass, such a per- 
turbed black hole in general relativity quickly settles back 
down to the stationary Kerr solution and in the process 
emits characteristic gravitational ringdown radiation in 
the form of so-called quasi-normal modes. These modes 
are uniquely determined by the mass and spin of the 
black hole and a detection of at least two of these modes, 
therefore, tests whether the end-state of the merger is 
a Kerr black hole constituting another strong-field test 
of general relativity (SJ) . Electromagnetic tests also in- 
clude the monitoring of close stellar orbits around Sgr A* 
135] and of pulsar black-hole binaries [36] , which are both 
weak-field tests. 

Strong-field tests of the no-hair theorem in the elec- 
tromagnetic spectrum are based on the imprints of the 
underlying spacetime on the emission from the accre- 



tion flows around black hole candidates, which can al- 
ter the observed signals in a characteristic manner [16|. 
Possible observables include the direct imaging of super- 
massive black holes with very-long baseline interferomet- 
ric techniques [ITij . relativistically broadened iron lines 
[H [13, nil, quasi-periodic variability E, IM HI, ther- 
mal accretion disk spectra 
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X-ray polarization 



[37j , and, potentially, jets 

In many of these tests, accretion flows are modeled as 
thin disks that lie in the equatorial plane of the compact 
object and that can be either optically thin or optically 
thick. In such models, the particles in the disk typically 
move on circular orbits. Deviations from the Kerr metric 
can have a strong effect primarily on the flux emitted 
in the part of such disks that is closest to the compact 
object, because the spacetime curvature experienced by 
the disk is largest at the inner edge. The inner edge 
of a disk around a black hole candidate is many times 
taken to coincide with the ISCO. Since the temperature 
of a thermal disk spectrum reaches its maximum near the 
ISCO (e.g., [3^), this inner region of an accretion disk 
can contribute a significant portion to the total observed 
disk flux. 

While the location of the ISCO in the Kerr metric 
(measured in units of the mass M) depends only on the 
value of the spin, the location of the ISCO usually de- 
pends on both the spin and the deviation parameters in 
a Kcrr-like metric (e.g., [lal)- At the ISCO radius, the 
radial epicyclic frequency vanishes, and circular equato- 
rial orbits are radially unstable inside of the ISCO (the 
so-called plunging region). In certain Kerr-like metrics 
such as the quasi-Kerr Q and Manko-Novikov [2(| met- 
rics, however, for certain values of the deviation parame- 
ters, all circular equatorial orbits are radially stable. In- 
stead, the vertical epicyclic frequency vanishes at a radius 
r ~ a few M , and such orbits experience an instability 
in the vertical direction [ll|, [l^ . Very close to the ori- 
gin, this vertical instability may di sapp ear, and circular 
equatorial orbits are, again, stable |11| . which may lead 
to characteristic gravitational-wave emission during an 
extreme mass-ratio inspiral [ij, |3^. As was shown in 
Ref. [40|, the vertical instability of circular equatorial 
orbits in the Manko-Novikov metric can cause geometri- 
cally thin accretion disks to become geometrically thick 
in the innermost region and the disk plasma may get 
trapped inside of the inner disk edge. The orbital sta- 
bility of these particles, in turn, depends sensitively on 
their energies and angular momenta, which may lead to 
an outflow of disk plasma in the direction vertical to the 
disk nil,!!!. 

In this paper, I study the properties of circular equato- 
rial orbits and the inner edges of thin accretion disks in 
the Kerr-like metric constructed by Johannsen & Psaltis 
|27| . This metric depends on a set of free parameters and 
derives from a general modification of the (t, t) and (r, r) 
components of the Schwarzschild metric with a successive 
transformation of the Newman- Janis type [43| . In general 
relativity, the Newman- Janis transformation is a com- 



plex rotation that transforms the Schwarzschild metric 
into the Kerr metric [4j) . While this transformation gen- 
erally does not relate static solutions to the correspond- 
ing stationary solutions in modified gravity theories, it 
can be used to introduce spin to modified Schwarzschild 
metrics (see, also, Refs. [Ij, [IJ). If all deviation pa- 
rameters vanish, such a transformed static metric, then, 
reduces smoothly to the Kerr metric. Contrary to other 
Kerr-like metrics such as the quasi-Kerr Q and Manko- 
Novikov [26[ metrics, which often contain singularities or 
closed timelike curves outside of the central object that 
can limit the utility of a given metric to model accre- 
tion flows around black hole candidates, the metric by 
Johannsen & Psaltis is free of such pathologies [23, |23| . 

In this metric, an ISCO exists for most deviations from 
the Kerr metric, while circular equatorial orbits become 
vertically unstable if an ISCO does not exist 27|. In the 
former case, the location of the ISCO as well as certain 
properties of circular equatorial orbits were already ana- 
lyzed in Ref. [23. Here I focus on the case when an ISCO 
does not exist. First I derive expressions of the energy 
and axial angular momentum of a particle on a circular 
equatorial orbit as well as the radial and vertical epicyclic 
frequencies. I calculate the radial locations, at which cir- 
cular equatorial orbits become unstable in the vertical 
direction, and show that these locations can be shifted 
significantly as a function of the spin and the deviation 
from the Kerr metric. I compare these radial locations 
to those of the ISCO in the part of the parameter space, 
where an ISCO exists. 

As a possible observable of accretion disks, whose inner 
edges are determined by the vertical instability of circu- 
lar equatorial orbits instead of the location of the ISCO, 
I simulate pro files of relativistically broadened iron lines 
(see Refs. [2j, [2^ for a discussion of iron lines in this 
metric emitted from disks that terminate at the ISCO). 
I show that these profiles can depend significantly on the 
deviation from the Kerr metric, which affects primarily 
the flux ratio of the two peaks in the double-horn shape of 
the line profile as well as the low-energy part of the line. 
While at a given value of the spin the ISCO radius de- 
creases for increasing values of the deviation parameter, 
the radius at which the vertical instability occurs likewise 
increases. This effect results in a different dependence of 
the profile shape on the deviation parameter and, in par- 
ticular, in different flux ratios of the two peaks in iron 
lines with unequal values of the deviation parameter for 
which the ISCO and the radius of the vertical instability 
are located at the same radius. I also show, however, that 
proflles for different choices of the spin and the deviation 
parameter that correspond to the same location of the 
inner egde as determined by the vertical instability are 
very similar. 

In Section |lTl I discuss some of the properties of the 
metric [27J and derive the energy and axial angular mo- 
mentum of a particle on a circular equatorial orbit. I 
analyze the stability of inner accretion disk edges in this 
metric in Section IIIIl In Section IIVI I simulate iron line 
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FIG. 1. Different shapes of the (singular) Killing horizon for values of the spin 



0.8M and deviation parameter 63. Left: 



£3 — 0.85 < e3°""'^(0.8A/), where e3°""''(a) denotes the boundary between the regions of the parameter space with different 
topologies of the Killing horizon. The inner and outer Killing horizons have spherical topology. Center: 63 = 0.92799 « 
^boundj-Q gjy^-j The inner and outer Kilhng horizons merge in the equatorial plane. Right: £3 = 1.0 > e3°""'*(Q.8A/). The inner 
and outer Killing horizons have split into two sphere-like surfaces located above and below the equatorial plane, respectively. 
For all values of the parameter £3, the origin is likewise singular. 



profiles, and I discuss the results in Section |Vl Through- 
out this paper, I use geometric units and set G — c = 1, 
where G and c are the gravitational constant and the 
speed of light, respectively. 



II. METRIC PROPERTIES 

In this section, I discuss some of the properties of the 
metric constructed by Johannsen & Psaltis [23] and cal- 
culate the energy and axial angular momentum of a par- 
ticle on a circular equatorial orbit. 

In Boyer-Lindquist-like coordinates, i.e., in spherical- 
like coordinates {t,r,9,(f>) which reduce to the standard 
Boyer-Lindquist coordinates if all deviations from the 
Kerr metric vanish, the metric is given by the expres- 
sions 



gu = ^[l + h{r,9)] 1- 



2Mr 



2aMrsm'^9, 

gt4, = ^^ [l + h{r,t 



9rr 



T,[l + h{r,e)] 



A + a^sm^0h{r,e)' 
gee = S, 



sm 



2 , 2 
r + a 



2a'^Mrsin^d 



h{r,e) 



a^{T. + 2Mr)sin^e' 



where 



A = r"' - 2Mr + a^ 
Y. = r^ + 0!^ cos^ 9, 



(1) 



(2) 
(3) 



and where 






fc=i 



(4) 



is a function that contains the deviation parameters efc. 
In this paper, I will use this metric with only one non-zero 
parameter, so that the function h(r^ 9) reduces to 



h{r,9)^e:, 



M'^r 



(5) 



In the limit £3 — > 0, this metric reduces to the Kerr met- 
ric. 

As shown in Ref. 28], if £3 7^ 0, this metric harbors a 
naked singularity, which is located at the Killing horizon. 
For small deviations from the Kerr metric, if the metric 
elements given by Eq. ([IJ are expanded to first order 
in the parameter £3, this singularity disappears, and the 
metric describes a black hole. 

While the Killing horizon has spherical topology for 
most values of the spin and the deviation parameter, for 
values of the parameter £3 > ground ^ where 



d/„\ — 



\a) 



1 



'3125(a/M)2 
-160(a/M)2 



1024 ^4 + v/16 - I5{a/M)A 



40 + 7v/l6 - 15(a/M)' 



+150(a/Af )" ( 15 + VI6 - 15(a/Af)- 



,(6) 



the topology of the Killing horizon is disjoint consisting of 
two spherical surfaces located above and below the equa- 
torial plane and centered at the symmetry axis (except 
for a = zfcM, which corresponds to the extremal Kerr 
black holes) [2^ . In Figure [1] I plot the location of the 



Killing horizon in the xz-plane, where a; = \^r^ + o? sin0 
and z = rcos9, for a value of the spin a ~ 0.8M and 
several values of the parameter 63 . As shown in Ref . [23] , 
the topology of the Killing horizon is also related to the 
existence of an ISCO in this spacetime. An ISCO exists 
for values of the spin a < as well as for values of the 
spin a > and values of the parameter £3 < eJ^°""d (as 
long as £3 < 30). 



The properties of a particle on a circular equatorial 
orbit depend directly on its orbital energy and angular 



momentum. For the case when £3 < e 



bound 



these ex- 



pressions were calculated in Ref. [27| and were written 
in terms of several functions Pi, i = 1,...,6, which are 
polynomials in the radial coordinate r apart from a few 
square root terms. As was pointed out in Ref. [37| , a sign 
change can occur in the axial angular momentum Lz at a 
radius very close to the central object, where the function 
P5 has a root. Circular equatorial orbits likewise exist in 
the part of the parameter space, where £3 > e^°^"<i. In 
the following, I generalize the expressions of the energy 
and axial angular momentum of a particle on a circular 
equatorial orbit that were derived in Ref. [23] to include 
the entire parameter space. The resulting expressions of 
the energy and axial angular momentum of particles on 
such orbits take the full root structure of the functions 
Pi to Pg into account. 

The derivation is very similar to the one in Ref. [27|. 



Denoting the 4-momentum of the particle by 

dx" 



P 



M- 



dT 



(7) 



where fi is its rest mass, the energy E = —pt and ax- 
ial angular momentum Lz = Prf, are conserved along its 
orbit. Solving the equation 



P Pa 



-M 



(8) 



in the equatorial plane for the radial and vertical mo- 
menta and inserting the constants of the motion, I obtain 
the effective potential 



K.ff(r) ^ \^i^ 



dr\ 

9rr 1 ^ j 



-^g''E^~2g'^ELz+g^^Ll + ^i^), 



(9) 



which governs the motion of particles in the equatorial 
plane. 

Particle motion is only allowed if Vcfi{r) > 0, which I 
can express in terms of the constants of the motion. At 
a given radius r and for a given axial angular momen- 
tum Lz, particle motion is only allowed for values of the 
energy E > i?min- From the equation 



Kff(r-) = 0, 
I obtain the minimum energy 



(10) 



E„ 



r2[r6 + a'^{r + 2M){r^ + esM^)] 



2aMr^Lz{r^ + €j,M^) 



+ ^Jr^{r^ + £3M3)[r4(r - 2M) -t- a?{r^ -f ezM^)]{r^Ll + [r^ + a?{r + 2M)(r^ + e3M^)]fi^} 



(11) 



Smaller values of the energy are excluded, because these 
correspond to negative values of the energy as observed 
in the inertial frame of a local observer. Note, however, 
that the minimum energy (as observed at radial infinity) 
can be negative. 

Solving the equation 



dVeajr) 
dr 







(12) 



for the axial angular momentum and combining the result 
with Eq. ([T0| . I derive the expressions 



(13) 
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1 

-^1^:6 1 


1— 


+ P2 
P3 


Lz 


1 


cr2\ 


IP^ 


A^ 


r^Pe 


IP3 



+a36aM{r^ + e^M^) 



P 



P3 



(14) 



I 

where cri_3 = ±1. These expressions are valid for all 
values of the spin |a| < M and the parameter £3, provided 
that a circular equatorial orbit exists. 

The signs of the constants cti_3 depend on the (real) 
root structure of the functions Pi to P5. In this paper, 
I only discuss the real roots of these functions if they lie 
in the range r > as well as outside of the ISCO, if an 
ISCO exist. Only in this case are these roots relevant 
for the discussion of circular equatorial orbits. It turns 
out that the roots of the polynomial Pg , which occurs in 
the denominator of the axial angular momentum, play 
no role. If such a root exists in the range that is relevant 
for circular equatorial orbits, the numerator of Eq. p^ 
likewise vanishes, and the expression of the axial angular 
momentum is smooth across the root. As an example of 
the roots of the functions Pi to P5, I plot the function 
P^ir) in Figure [5] for values of the spin a = 0.7 M and 
the parameter £3 = 2.0. The function depends very sen- 
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FIG. 2. Function P5{r) as a function of the radius r with 
parameters a — 0.7M and es = 2.0. Top panel: The function 
Ps (r) is very steep and has several local minima and maxima 
as well as five roots marked by red vertical dashed lines, one 
single root labeled "ns" and two double roots labeled "nsa" 
and "nsb", respectively. Lower panel: Enlarged view of the 
function Ps around the roots on the r-axis. 



sitively on the radius r and is very steep at most radii. In 
addition to several local maxima and minima, the func- 
tion contains five roots, four of which are double roots. 

For most circular orbits, the expressions of the energy 
and axial angular momentum given by Eqs. (|13p and (J14p 
are identical to the ones in Ref. [23| (apart from a minor 
rearrangement of the function P3), i.e., cri = +1, (72 = 
± 1 , and 0-3 = =F 1 , where the upper signs refer to prograde 
orbits, while the lower signs refer to retrograde orbits. 
For values of the parameter —10 < £3 < 10, sign changes 
can only occur at radii r ^ 4.5M, most of which lie 
inside of the ISCO if an ISCO exists. If an ISCO exists, 
these sign changes usually only affect the energy and axial 
angular momentum slightly. 

If eg < eb°""d^ cTi = +1. If eg > ebound^ ^l^g yj^l^g 

of fJi depends on the existence of a root of the function 
Pi + P2- The values of a2 and 0-3 depend on the root 
structure of all functions Pi to P5 , and several different 
cases have to be considered. Outside of the Killing hori- 
zon, if the function P3 has roots, the functions Pi + P2 
and P5 have roots that coincide with the roots of the 
function P3, and the expressions of the energy and ax- 
ial angular momentum remain finite at the roots. If the 
Killing horizon has spherical topology, the expressions 
of the energy and axial angular momentum become in- 
finitely large at the singularity. If the Killing horizon is 
disjoint and there is no naked singularity in the equa- 
torial plane apart from the origin, the energy and axial 
angular momentum change sign at a radius r '^ M , and 
E, Lz —T' —00 as r — )■ 0. In Figure[31 1 plot the energy and 




FIG. 3. (Top) Energy E and (bottom) axial angular mo- 
mentum Lz in units of the rest mass of a particle on a 
circular equatorial orbit as observed at radial infinity for 
a spin value a = 0.8A/ and several values of the parame- 
ter £3. If the Killing horizon has spherical topology [i.e., if 
£3 < e3°""'*(0.8M) ^ 0.92799], both E and L^ diverge to plus 
infinity at the horizon and are imaginary inside. Otherwise, 
they change sign near the outermost projected location of the 
Killing horizon onto the equatorial plane and diverge to minus 
infinity at the origin. 



axial angular momentum for a value of the spin a — 0.8M 
for several values of the parameter €3 , where the Killing 
horizon has spherical topology for values of the parame- 
ter £3 < ground ^ 0.92799 and disjoint topology for values 
of the parameter £3 > £1^°""^^ j jjg^ ^-j^^^ ^^y^ expressions of 
the energy and axial angular momentum given the roots 
of the functions Pi to P5 in Appendix A. 

Note that in certain cases both the energy and axial 
angular momentum as observed at radial infinity can be 
negative. While negative values of the axial angular mo- 



mentum simply correspond to retrograde orbits, negative 
values of the energy are best interpreted in the inertial 
frame of a local observer (see discussion in Ref. [4J|). 
For such a frame I can choose an orthonormal basis with 
vectors 



negative values of the energy can only occur if 

^2 



1- 



5t0 - 9tt9^^ 



>0, 



(23) 



which obey the relations 

6 a^ /jdp.i' — VaP^ 



(15) 



(16) 



where 77^3 = diag( — 1, 1, 1, 1) is the Minkowski metric. 
Making the ansatz for the basis vectors 




(17) 



I use the relations given by Eq. ([T6| to solve for the pa- 
rameters 7] and 7. I obtain the solutions 



fft0 - 91194,4, 



I r2[r6 + a2(r + 2M)(r3 + e3M3)] 

(r3 + e3Af3)[r4(r - 2Af) + a^ir^ + esM^)] ' 



94>4> 



94>4> V 9t4> - 9tt94>4> 



'2aMr 



esAP 



[r6 + a2(r + 2A/)(r3 + e^A'P)] 



[r'^ir - 2M) + a^ir^ + esM^)] 



(18) 



The locally measured energy and axial angular mo- 
mentum are the projections of the 4-momentum of the 
particle onto the basis vectors e^ and e t given by Eq. ([TTl) : 



(19) 
(20) 



Solving these equations for the energy, I can write the 
energy as 




(21) 



Since all physical particles must have momenta with 



P'> P"" , 



(22) 



9tt > 0, 



(24) 



where I used Eqs. p8|) . Therefore, negative values of the 
energy as observed at radial infinity can only occur in the 
ergosphere. 

In Figure m plot the energy as a function of the radius 
for values of the spin a = 0.8Af and the parameter £3 = 
0.92799 as observed in the local frame given by Eqs. (|17l) . 
While the energy observed at radial infinity is negative 
at small radii (c.f.. Figure [3|), it is always positive in the 
local frame. 
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FIG. 4. Energy iJiocai in units of the rest mass of a particle on 
a circular equatorial orbit as observed by a local observer for 
values of the spin a = 0.8Af and the parameter 63 = 0.92799. 
The energy is positive at all radii r > 0. In particular, at 
least at the radii shown in this figure, the energy is greater 
than or equal to the rest mass of the particle, where equality 
holds only at a radius r ~ 1.31Af (see inset). 



III. INNER ACCRETION DISK EDGES 



For values of the spin a < as well as of the spin a > 
and the parameter £3 < e^ound^ circular equatorial orbits 
are radially unstable inside of the IS CO radius, while cir- 
cular equatorial orbits are radially stable for values of the 
spin a > and the parameter £3 > ground j^ ^^^ latter 
region, circular equatorial orbits are, instead, unstable 
against small perturbations in the direction vertical to 
the equatorial plane [23|. At this radius, accretion disks 



are likely to acquire a certain thickness and may be dis- 
rupted (see Ref. pOl - t42| ). Near the origin, the vertical 
instability disappears, and a second region of stable cir- 
cular orbits emerges inside of this radius, which extends 
down to the origin. In principle, this region could harbor 
a small inner part of the accretion disk, which is discon- 
nected from the accretion disk located at radii outside of 
the vertical instability (see Refs. jjU, 39, 45]). 

In order to determine the radial locations, where cir- 
cular equatorial motion is either stable or unstable in 
the radial or vertical directions, I calculate the radial 
and vertical epicyclic frequencies fir and fie following 
the derivation in Ref. [18|]- From the effective potential 
given by Eq. pHll . radial and vertical motions around a 
circular equatorial orbit at a radius tq are governed by 
the equations 



d9_V 
'di) 



Kff 



ffrr(p*)^ 



V. 



ff' 



v: 



ffj 



(25) 
(26) 



respectively. 

I, then, introduce small perturbations &r and bd and 
take the coordinate time derivative of Eqs. (1251) and (l26t . 
which yields 






eff 



dr'^ 



Sr, 



d^iS9) _ d^V^ 



cff. 



dt^ 



d9^ 



(27) 
(28) 



From these expressions, I derive the radial and vertical 
epicyclic frequencies as 






dr^ 
d0^ 



(29) 
(30) 



where the second derivatives are evaluated at r — vq. 
These expressions are lengthy, and I do not write them 
here explicitly. 

In Figure [5l I plot the radial and vertical epicyclic fre- 
quencies versus the radius for a spin value a = 0.8M 
and several values of the parameter e^. For values of the 
parameter £3 < 6^°^""^ ^ 0.92799, circular orbits are ra- 
dially unstable inside of the ISCO radius and are stable 
at and outside of the ISCO radius. For values of the pa- 
rameter £3 > gbound^ ^Y\ circular orbits are radially stable 
but are vertically unstable over a small range of radii, 
where the vertical epicyclic frequency is imaginary. 

In Figure [SI I plot contours of constant ISCO radius 
as a function of the spin and the deviation parameter. 
An ISCO exists in the part of the parameter space where 
a < and where £3 < gbound^ a > 0. Contours of con- 
stant ISCO radius decrease for increasing values of the 
spin and decreasing values of the parameter £3 [271 ] . In the 
part of the parameter space, where £3 > e^ound^ a > 0, 
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FIG. 5. Radial (top) and vertical (bottom) epicyclic frequen- 
cies Ur = ilr/2'K and ve = fle/2TT for a lOM0-black hole for 
a — O.SM and several values of the parameter £3. If an ISCO 
exists, the radial epicyclic frequency vanishes at the ISCO ra- 
dius, and circular equatorial orbits are radially unstable inside 
of this radius. If an ISCO does not exist, all circular orbits are 
radially stable but vertically unstable roughly in the interval 
O.IM <r< 1.7M. 



I plot contours of the (outer) radii, at which the verti- 
cal epicyclic frequency vanishes causing circular equato- 
rial orbits to become unstable in the vertical direction. 
These contours increase for increasing values of the spin 
and increasing values of the parameter £3 . The transition 
between radial and vertical instability of circular orbits 
across the boundary marked by the function £3°""'^ (a), 
a > 0, is generally not continuous, i.e., these instabil- 
ities occur at different radii just above and below this 
boundary. 

In the region of the parameter space, where an ISCO 
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FIG. 6. Contours of (blue curves) constant ISCO radius and 
(green curves) constant radius, at which the vertical epicyclic 
frequency vanishes, if no ISCO exists. The boundary be- 
tween these regions (red curve) is determined by the function 
(a), a > 0, given by Eq. ([6|. 



bound 
•^3 



exists, circular equatorial motion is also unstable in the 
vertical direction. This instability, however, usually oc- 
curs inside of the ISCO and does not affect the struc- 
ture of accretion disks. In a small and narrow region 
below the boundary e3°™'^(a), a > 0, however, the ra- 
dius, where circular motion becomes vertically unsta- 
ble, lies slightly outside of the radius, where circular 
motion becomes radially unstable. This region lies in 



the spin range 0.57M 



< 



< M and has an extent 



Ae3(a) 



— ^bc 



Ka) - e3,isco(a) < 0.03, where e3,isco(a) 



is the value of the parameter £3 at which the locations of 
the ISCO and the vertical instability coincide at a given 
value of the spin. I plot this region in Figure [71 In Ta- 
ble I, I summarize the properties of the inner edges of 
accretion disks and the different topologies of the Killing 
horizon. 

Killing Horizon and Inner Accretion Disk Edges 
Spin Deviation Killing Horizon ISCO Vertical 

Parameter Topology Instability 



a < £3 < eoounti spherical 
3 disjoint 



■^3 

a < 0^ £3 > fbound 



a > €3 < fbound spherical 

g > 0" £3 > gbound 



disjoint 



yes inside of ISCO 

yes inside of ISCO 

yes inside of ISCO'" 

no yes 



'^ Kerr black hole if a = ±M with a spherical Killing horizon and 
the ISCO located at r = M and r = 9M, respectively. 

^ Except for values of the spin and deviation parameter in a 
narrow region along the boundary £3°'""^ (see Figure [TJ. 



TABLE I. 
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FIG. 7. Top: Regions, where an ISCO exists (labeled "ra- 
dial instability") and where all circular equatorial orbits are 
radially stable but vertically unstable (labeled "vertical insta- 
bility"). The boundary between these regions (blue curve) is 
determined by the function e3°"°'^(o), a > 0, given by Eq. ((6)). 
If an ISCO exists, circular equatorial motion is also vertically 
unstable, but this vertical instability only occurs inside of 
the ISCO except for values of the spin and the parameter 
£3, which lie in a very narrow (red shaded) region along the 
boundary in the range 0.57M ^ o < M. Bottom: Thickness 
Ae3 of the red shaded region in the top panel. 



IV. RELATIVISTICALLY BROADENED IRON 
LINES 

As a potential observational signature of the radiation 
emitted from accretion disks whose inner edges are de- 
termined by the vertical instability, I simulate profiles of 
relativisitcally broadened iron lines. I use the algorithm 
already developed in Ref. [2^ with the proper implemen- 
tation of the energy and axial angular momentum in the 



range where £3 > £b°""'i as listed in Appendix A. 

As in Ref. [23i], I place an observer at a large distance 
d from the central object and at an inclination i from 
the rotation axis. I assume that the accretion disk is 
geometrically thin and optically thick and that it has a 
radial extent from the (outer) radius, where the vertical 
epicyclic frequency vanishes, to some outer radius Tout- I 
assume that each particle in the disk moves on a circular 
Keplerian orbit, i.e., with the local Keplerian velocity 
u = (w*,0,0,w'^), where [H 



9ttg4,4, 



-9tci> 
9ttLz 



9t^ 



(31) 



(32) 



9tt9<P<P 

Further I assume that the emission from the disk is 
monochromatic in the restframe with a typical energy 
Eq w 6.4 keV and has an isotropic emissivity e(r) ex r~", 
where a is the emissivity index. 

I calculate the photon paths to the distant observer via 
a fourth-order Runge-Kutta integration. These photons 
experience a redshift 
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Eim 






(33) 



where the subscripts "im" and "d" refer to the image 
plane and the accretion disk, respectively, and where fc'* 
is the 4- momentum of a given photon. The observed total 
specific flux is then given by the expression |23| 



dao / d/3o/(ao,/3o)(5[£'p - £;og(ao,/3o)], 



(34) 

where ag and Pq are the (Cartesian) coordinates of the 
image plane and where I{ao,l3o) and Ep are the specific 
intensity and energy, respectively. 

In order to avoid any contact with the naked singu- 
larity located at the Killing horizon, I excise all radii 
r < 1.05 X Thor, where rhor is the maximal radius of the 
Killing horizon. This cutoff acts as an artificial event 
horizon, and any particle passing through this horizon 
is considered "captured" and does not contribute to the 
observed line flux. In principle, the flux emitted from the 
potential inner part of the accretion disk around the ori- 
gin could also be included. In practice, however, if such 
an inner disk exists, the electromagnetic radiation emit- 
ted from this inner region likely accounts for only a small 
fraction of the total flux and and is highly redshifted. 

In Figure ISl I plot iron line profiles for values of the spin 
a = 0.9M, index a = 2, rout = 15M, i = 30°, and several 
values of the deviation parameter. These lines have the 
characteristic double-horn profile with a "blue" peak at 
the energy E Ri Eq and a "red" peak at a slightly lower 
energy as well as the "red tail" at low energies which 
terminates at the energy with the highest gravitational 
redshift that corresponds to the emission from the inner 
accretion disk edge. The effect of changing the parame- 
ter €3 is similar to its effect analyzed in Ref. [27], which 
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FIG. 8. Profiles of relativistically broadened iron lines in units 
of the rest frame energy Eo for values of the spin a — 0.9AI, 
emissivity index a — 2, outer disk radius rout = 15i\f , disk 
inclination i — 30° , and several values of the deviation param- 
eter. In the case of a Kerr black hole (ea = 0), the inner edge 
is located at the ISCO. For values of the parameter 63 = 1.0 
and £3 = 5.0 the inner edge is located at the radius, at which 
the vertical epicyclic frequency vanishes. Changing the pa- 
rameter £3 primarily affects the flux ratio of the two peaks of 
a given profile as well as the low-energy spectrum correspond- 
ing to the location of the inner disk edge. 



focused on values of the deviation parameter £3 < ground ^ 
For increasing values of the parameter £3, the flux of the 
blue peak diminishes, while the flux of the red peak in- 
creases and can even be larger than the flux of the blue 
peak, and the end of the red tail shifts according to the 
shift of the inner disk edge. In particular, changing the 
parameter £3 affects the flux ratio of the two peaks, which 
is a potential observable. 

At a given value of the spin, the flux ratio of the two 
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peaks can be significantly different for iron lines with 
different values of the deviation parameter, for which the 
ISCO and the radius of the vertical instability lie at the 
same location. In Figure [HI I plot iron lines with a spin 
parameter a — 0.9M and values of the deviation param- 
eter 63 = —0.0775 and 63 ~ 5.0, for which the inner 
accretion disk edge as determined by either the ISCO 
and the vertical instability, respectively, is located at a 
radius r w 2.423M. The flux ratio of the peaks is much 
larger in the case if the inner edge is set by the ISCO 
(£3 = —0.0775) than in the case if the inner edge is set 
by the radius of the vertical instability (es — 5.0). In 
addition, in the former case the red tail is much more 
elongated, while it falls off abruptly in the latter case. 




X 2 



- 


a = 0.9M, a = 2, r„„j=15M 




- 


£3=5.0 i\y 

e3=-0.0775 r^ 


- 


- 


1 1 J--^ / 1 1 1 1 1 1 


1 1 



0.5 



X 2 
3 







_ a = 2. r,„^=15M 


- 


. a = 0.60M. £3=4.892 ft 


_ 


- a = 0.99M, €3=2.0 £ / 


- 


i W^ 


_ 


/ 


- 


-J^ , 1 





0.5 



E/Eo 



E/E„ 



FIG. 9. Profiles of relativistically broadened iron lines in units 
of the rest frame energy Eq for values of the spin a — 0.9M, 
emissivity index a — 2, outer disk radius rout = 15A/, disk in- 
clination i — 30° , and two values of the deviation parameter. 
In both cases, the value of the deviation parameter is chosen 
so that the inner edge of the accretion disk as determined 
by either the ISCO (ea = —0.0775) or the vertical instabil- 
ity (q = 5.0) is located at the same radius (r « 2.423M). 
The flux ratio of the peaks is much larger in the former case 
than in the latter case and the red tail extends to much lower 
energies. 



As shown in Ref. [23], profiles for sets of parameters 
with the same location of the ISCO are practically indis- 
tinguishable. In Figure [TUl I plot line profiles with two 
sets of parameters that correspond to the same radius of 
the vertical instability. In this case, the profiles are like- 
wise very similar indicating a strong correlation between 
the locations of the inner edge and the line profiles. 



FIG. 10. Iron line profiles in units of the rest frame energy Eq 
with values of the spin and the deviation parameter chosen so 
that the inner disk edges coincide. The line profiles are very 
similar. 



V. DISCUSSION 

In this paper, I calculated the energy and axial angular 
momentum of a particle on a circular equatorial orbits as 
well as the radial and vertical epicyclic frequencies in the 
Kerr-like metric constructed by Johannsen & Psaltis 27]. 
These functions depend explicitly on a set of functions Pi 
to Pq and their root structure, which I analyzed in detail. 

I investigated the properties of the inner edges of thin 
accretion disks in this spacetime. While the ISCO is lo- 
cated at the radius, where the radial epicyclic frequency 
vanishes, an ISCO does not exist and circular equato- 
rial orbits are radially stable at all radii in the range 
£3 ^ gbound^ where this lower bound on the parameter e^ 
is given by Eq. ^. Instead, circular equatorial orbits be- 
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come unstable in the vertical direction. I calculated the 
radii at which this instability occurs and showed that 
these locations depend significantly on the spin and the 
deviation parameter. The vertical instability exists even 
in the case if an ISCO does exists, but I argued that 
the vertical instability plays no role in the structure of 
thin accretion disks in this case, because it always lies 
inside of the ISCO unless the deviation parameter £3 lies 
in a narrow region along the boundary £3°^"'^ (a), which 
I calculated explicitly. 

As a possible observable of accretion disks, which ter- 
minate at the radius where the vertical instability occurs, 
I simulated profiles of relativistically broadened iron lines 
and showed that these profiles depend significantly on the 
values of the spin and the deviation parameter. In par- 
ticular, I showed that the flux ratio of the red and blue 
peaks can be significantly altered, which can serve as a 
probe of the deviation from the Kerr metric and which 
distinguishes iron lines with the same spin but unequal 

I 



values of the deviation parameter for which the inner disk 
radius as determined by the ISCO and the vertical insta- 
bility, respectively, is identical. As in the case of disks 
that terminate at the ISCO (see Ref. 23]), I also showed, 
however, that profiles for different choices of the spin and 
the deviation parameter that correspond to the same lo- 
cation of the inner egde as determined by the vertical 
instability are very similar. 
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Appendix A: Energy and Axial Angular Momentum 

The functions Pi, i = 1, ..., 6, that occur in the expressions of the energy and axial angular momentum of a particle 
on a circular equatorial orbit given by Eqs. (fT3| and (fT4|) have the form 

-r^ [2e3M^r^ (Sr^ - 8M^) + eJM^ {40M^ - 48Afr + 15r^) + 4r''(3r - 5M)] I, 
P2 ^21 2r'^ (r^o =f AIP4) + Mr^^ I 2r^ (-121/^ + \QMr - 7r^) 

+e3M^(r - 2Mf [e^Af^(5r - 12M) - Qe:iM^r^{hM - 2r) - 3r'^(8M - 3r)] I I, 

P3 =r* (l2e3M'' - be^M^r + QMr^ - 2r*f - Sa^Af (53 A/^ + r^f {besM^ + 2r^) , 
J —aB if ^41 < r < ^42 

Pi — ) J 

I aB else 

P5 ^M (r^ + esM^)^ 1 12e3a^Af^ {e3M^ - 2r^)^ {esM^ + r^)^ + a^r'' (esM^ + r^)^ 

{-AOelM^^ + AOeJM^'^r - ISe^^A/^r^ + 128e^A.f ^^r^ - 296e^A/^r* + die^M^r^ - 924ep//V^ 
+276elM^r'^ - 36e^AfV^ - SSOeaM^r^ + SOAesM^r^^ - 24e3Af2r" - 112Mr^^ + 16r^^) 
-2aV^ [48e^Af 1^ - 12e^Af ^^r - 52elM^^r^ + Sej^Af ^^^(Ses + 88) - 720e;^Af "r"* + 298e^A/"iV^ 



,2,^7^10 



(Al) 

(A2) 
(A3) 
(A4) 



-3e^Af ^V*'(13e3 + 480) + 516e;^Af ^'V^ + 2eiArr'^ - 6e^AfV^(3e3 -I- 508) -I- 2292e^Af 

-628e^AfV" -I- UesM^r^^ibes - 134) + llSSesM^r^^ - 296e3M^r^'^ + 2AM'^r^^{ej, - 9) + 120Afri^ 

-16r"] - r^^ {e^M^ + 6Mr^ - 2r^Y 

(96e^Af^ - IQeJM'^r + 15epf V^ + 72e3A/^r3 - UesM^r"^ + &e:>,M^r^ + 12A/r^ - 4r^) l 



-4P4 



(63 A/3 - 2r^f {esM^ + r^) + Qe^M^r^ {esM^ + 6Mr^ - 2r 



f4„2 



r3„3 



Pg ^^eiAr - eesM^r^ + esM-^r-" - dMr" + 2r^ 



(A5) 
(A6) 
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where 



B 



^M {e^AP + r^f {Qela^A'P + IGesMV^ - eegMV^ + 4r7) [a^ {e^AP + r^) + r^[i 



2M)X 



(A7) 



and where 7141 and 7142 denote the smaller and greater root of the function P4, respectively, if real roots exist for 
r > 0. Note the slight redefinition of the functions P4 and P5 with respect to the one in Ref. 27 1. 
Defining 



E+ = — 



L, 



1 P1+P2 



PPe 



P3 

fPr. 



6aAI{r^ + egAf ^) 



P1+P2 
P3 



L_ = -- 



1 



"3 



esM') 



P1+P2 



P3 



I can write the energy and axial angular momentum in the following manner: 

E+ if r > ni2 



E = 



if r < n^ 



E, 



—E+ else 



II £3 < £3 



(A8) 
(A9) 

(AlO) 
(All) 



L_l- if r > max(n3 , 715 ) 



if r > ^5 < 



else < 



if n-i = nz< 



A if r < n^a 
iir<n:i{ , \L+ r > ni2 



else 



else ■ 

I — L_ else 

L_ if sgn+ := sgn~ 



else < 





(A12) 



if r < n5a , ^ 

L+ else 

^ L^ else 

i+ if sgn+ = sgn^ 

else { I L 1 if 7^3 = n5 

else < 



where 



A = 



sgn"*" = sgn 



if 77> 



— i_ else 



"12 




(A13) 
(AM) 



, Li else 



In these expressions, ni2 denotes the root of the function Pi + P2, 713 is the largest root of P3, and 775, n^a, and 
7i5b denote the roots of P5, where n^ > n^a > n^b- Further, max(773,7i5) denotes the maximum of the roots n^ and 
775, and sgn^ denote the signs of the function P5 evaluated in a vicinity of the root 715 at radii r = 715 + Ar and 
r = 775 — Ar, respectively, where Ar <^ AI. The latter definition is a simple check for the existence of either a single 
or a double root at 715, which can easily be implemented in a numerical algorithm without having to compute the 
derivative of the function P5 . These expressions of the energy and axial angular momentum are valid for all values of 
the spin \a\ < AI and the parameter £3 if a circular equatorial orbit exists. 
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